Flow characterization inside airways with unsteady breathing patterns

Flow through human airways is characterized by unsteady flows, with flow separations at airway bifurcations. The oscillatory nature of airflow, unequal durations of inhalation time (IT), and exhalation time (ET) can facilitate gas exchange in higher generations of the human airway. Normal respiratory rate (RR) in adults ranges between 10-15 breaths per minute (bpm). RR varies in exercise conditions, mechanical ventilation strategies such as high-frequency oscillatory ventilation (HFOV), metabolic activities and pathological state to facilitate alveolar gas exchange. Previous studies characterized flows at steady inhalation and exhalation through airways. The individual effects of varying inhalation duration and breathing flow rate on flow through airways remains unknown. Our study focuses on various unsteady breathing patterns inside idealized airway models. The goal of this study is to characterize the effects of unsteady internal airflow through idealized airway geometries. Various scenarios of unsteady breathing patterns were simulated in ANSYS software (ANSYS, Inc., Canonsburg, PA, USA) to characterize the fluid dynamics involved in such an unsteady airflow mechanism. The first study includes unsteady breathing patterns such as normal, moderate, and high-frequency ventilation were investigated with variation in inhalation time (IT) to breathing time (BT) ratio. The second study includes abnormal breathing patterns such as tachypnea (~ 1.5x increase in RR), bradypnea (~ 1.5x decrease in RR), hyperpnea (deep breathing with abnormally large peak flow rate), and hypopnea (shallow breathing with abnormally low peak flow rate); and final study includes single nostril inhalation as in yoga pranayama breathing techniques. Simulations were performed for each breathing pattern as in internal airflow studies. Our results showed that secondary flow was an effective transport mechanism for flow inside idealized human airways. Airway local geometry plays a key role in flow distribution in higher generations. Discrepancy in the oscillatory flow relation Re/Wo^2 = 2L/D (L = stroke length; D = trachea diameter) was observed for IT/BT does not equal 50%, as L changed with IT/BT. We developed a modified dimensionless stroke length term including IT/BT. While viscous forces and convective acceleration were dominant for lower Wo, unsteady acceleration was dominant for higher Wo. Time to peak jet length during inhalation increased with an increase in breathing time. Single nostril and double nostril inhalation showed equal ventilation at higher generations in an idealized airway geometry